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Abstract 

Spacecraft observations of Saturn's rings show evidence of an active aggregation-disaggregation process 
triggered by periodic influences from the nearby moons. This leads to clumping and break-up of the 
ring particles at time-scales of the order of a few hours. A mathematical model has been developed 
to explain these dynamics in the Saturn's F-ring and B-ring [3], the implications of which are in close 
agreement with the empirical results. In this paper, we conduct a rigorous analysis of the proposed 
forced dynamical system for a class of continuous, periodic and zero-mean forcing functions that model 
the ring perturbations caused by the moon flybys. In specific, we derive the existence of at least one 
periodic solution to the dynamic system with the period equal to the forcing period of the moon. Further, 
conditions for the uniqueness and stability of the solution and bounds for the amplitudes of the periodic 
solution are derived. 

Keywords: Saturn ring, aggregation-disaggregation process, periodic solution, stability 
AMS2000 Subject Classification: 34C25, 85A99. 

1 Introduction 

In [3], Esposito et. al. have proposed a dynamical model to explain the aggregation and disaggregation 
processes observed in Saturn's F ring and the B ring outer edge due to the perturbation caused by the 
Saturn's moon, Prometheus and Mimas, respectively. Their numerical simulations demonstrate that the 
dynamical model explains certain spacecraft measurements corresponding to these regions in the Saturn's 
ring. Further, these numerical results strongly suggest the existence of a limit cycle for realistic values of 
the parameters. In this paper, we provide a rigorous mathematical proof of this fact as well as study other 
features that are relevant from a dynamical point of view. 

The model under consideration relates the mean aggregate mass M and the velocity dispersion v^^^ of 
ring particles at the above mentioned locations in the Saturn's ring, as follows 



dM _ M M 
~dt ~tZ.~\^J tZ. 
dvl vL{l-e^) M^ vjAMo) 
dt " T.„„ T.,., ^ ' 
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In the first equation of (ll.ip . the first term refers to the coagulation of the ring particles leading to a 
growth in the mean aggregate mass at an accretion rate (Tlcc is the accretion period) and the second 
term refers to the erosion (or break-up) of the larger ring particles due to collisions with other ring particles 
at coUisional rates T^on is the collisional period) when their velocity dispersion exceed Uth, the threshold 
velocity for sticking. 

In the second equation of (jl.ll) , the first term refers to the dissipative collisions between the ring particles 
leading to the reduction in velocity dispersion with e as the normal coefficient of restitution. The second 
term refers to the increase in the velocity dispersion to the escape velocity from mass M {vesc (Mq) is the 
escape velocity from an aggregate of mass Mq) due to the viscous stirring at a rate {Tstir is the stirring 
period), caused by the passage of the ring particles of mass M. Finally, the last term refers to the periodic 

forcing by the Saturn's moon causing disturbances in the velocity dispersion, where f{t) — ^^^"'^ cos(y^) 
describes the moon perturbation, where Tgyn is the forcing period, /? is the forcing amplitude. 

For convenience, we reduce the system to an equivalent dimensionless system with the following substi- 
tutions: X = -jj-^ y = and t = 77^, where T^i-h is the orbital period of the ring particles around Saturn, 

V^^^ J orb 

and pass to a more mathematical notation by using ^ . Then, the system in becomes 

x' = ax — bxy 

y' ^-cy + dx^-fit) (1.2) 

where 

^orb , Torh 7orb(l — £^) , Torbwic /^ n\ 

a=7f^,b^—,c= ,d=——^ (1.3) 

-'-acc J- coll -'-coll -'-stirL'th 

are all positive parameters, f (t) — Ao cos (ujt), Aq = w = ^j"^"'' = and T is the (normalized) 

forcing period. 

The above dynamical system has two main differences when compared to the predator-prey system. 
Firstly, in the second equation, we have the x'^ term instead of the coupling xy. Secondly, the external 
forcing is not standard in ecological models, since the typical forcing is parametrical (see for instance the 
model studied in |2]). 

Despite these differences, the aggregation-disaggregation model represented by this system can be well 
explained by drawing an analogy with the predator-prey system. Let the mean aggregate mass of the 
ring particles correspond to the prey population, and the velocity dispersion correspond to the predator 
population. The velocity dispersion 'feeds' off the accelerations from the aggregates' gravity (in reference to 
the second term of the second equation). As the velocity dispersion grows too large, it limits the 'prey' as 
high velocity collisions lead to to fragmentation of the ring particles (in reference to the second term of the 
first equation). In the absence of the interaction between the 'predator' and the 'prey', the 'prey' population 
(mean aggregate mass) grows and the predator population decays (in reference to the first terms of the 
two equations, respectively). In addition, the predator population is periodically controlled in a certain 
deterministic manner, as represented by the sinusoidal forcing function. 

In the following section, we conduct a thorough qualitative analysis of the above system and comment 
about the existence and nature of its solution, conditions on uniqueness and asymptotic stability of the 
solution and finally derive explicit bounds for the solutions. Along with these results, we derive useful 
insights about the Saturn's ring dynamics based on these results. 

Further, the results in this paper hold for any general forcing functions f{t) as long as it is a continuous 
periodic function of minimal period T = with zero mean value, i.e., Jp f{t)dt = 0. All the results will be 
presented for this general forcing function, and the sinusoidal forcing function will be explicitly studied in 
Section [51 



2 Existence of a periodic solution 

Numerical results in [3] strongly suggest the existence of a limit cycle for realistic values of the parameters. 
In the following theorem, we provide an analytical proof for the same. Further, the result holds for all 
feasible values of the system parameters. 
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Theorem 1 System (jl.2p has at least one T-periodic solution {x,y) with x{t) > for all t. 
Proof. By introducing the change z = Inx, system (jl.2p is transformed into 

z' = a — by 

(2.4) 

2/' = -c2/ + dcxp(2z)-/(t). 

The theorem is proved by showing that the above system is equivalent to a second order differential 
equation of the Duffing type, and then invoking the Landesman-Lazer conditions for the existence of a 
periodic solution. See Appendix \K\ for the complete proof. □ 

The above result leads to a physical interpretation that the periodic forcing by the moon gives a limit 
cycle in the mean aggregate mass and velocity dispersion with the same period as the forcing function. 

It is interesting to remark that the latter result can be extended to a family of forcing terms with more 
complicated recurrence including quasi periodic forcing. The main result in pi implies the existence of a 
quasi periodic solution and the boundedness of all solutions in the future. 

3 Uniqueness and asymptotic stability 

In this section, we derive conditions for the existence of a unique periodic solution and comment about the 
stability of the solution. 

Theorem 2 When 

acT<2+—, (3.5) 
system (jl.2p has a unique T-periodic solution, which is asymptotically stable. 

Proof. See Appendix [Bl □ 



Remark 1 Observe that for / = 0, the unique equilibrium of system (II. 2p is hyperbolic, hence robust to 
small perturbations. Therefore, if Aq is small enough, (|1.2p presents an asymptotically stable T-periodic 
orbit. The remarkable feature of the latter result is that stability is proved on a certain region of the involved 
parameters independently of the size of the external force f. 

As we will see in Section [5l Theorem [2] is applicable to the Saturn's F-ring, in fact for the F-ring 
ac ~ 16t^(1 — e^) ~ 0.104 so condition p.Sp is satisfied. The condition p.Sp is also satisfied for the Saturn's 
B-ring outer edge when the optical depth is r = .5. In the following theorem, a condition for the existence 
and stability of the solution is derived that depends on the nature of the forcing function via its supremum 
norm, denoted by II/Hq^. 

Theorem 3 Under the assumption 

2 2 

2ac + 2b\\f\\^<j + ^, (3.6) 
system (|1.2p has at least one asymptotically stable T-periodic solution. 

Proof. The above result is obtained by using the technique of lower and upper solutions [1]. See Appendix 
[Cl for the complete proof. □ 

In the context of the Saturn's ring dynamics, Theorem[3]and Theorem [3] show that the forced ring system 
has a limit cycle behavior around the fixed point and when the external forcing satisfies the above condition, 
it drives a stable limit cycle. Further, the above theorems provide a theoretical justification for the concentric 
limit cycles that were observed through numerical simulations in (3^ Fig. 16 and Fig. 17]. 
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4 Explicit bounds 



Theorem [T] provides a general existence result, but more concrete quantitative information about the location 
of the periodic orbits is desirable. In this section, we derive explicit estimates for the periodic solutions. In 
the next result, = max{/, 0}, /~ = max{— /, 0} denotes the positive and negative part respectively of a 
given function /, and ||/||p denotes the usual L^-norm. 



Theorem 4 Let {x,y) with x{t) > for all t a T -periodic solution of system (11.21) . Then, the following 
hounds hold 

v/fexp (-^^ 11/11,) < x{t) < yfexp (^^ ||/||,^ 
f (1 - cT) - < y{t) < f (1 + cT) + \\f-\\, 

for all t. 



(4.7) 



Proof. The bounds derived above are obtained by working with the second order differential equation of 
Duffing type (|1.12l) that is equivalent to the system in (|1.2I) . See Appendix |D] for the complete proof. □ 

Since the lower bound for x {t) is non-negative, the T-periodic solution x (t) is also non-negative for all 
t. The above argument acts as an alternate proof for the later part of the statement in Theorem [TJ Further, 
the bounds derived for the T-periodic solution for y (t) may not necessarily be tight. Especially, for small 
forcing amplitudes, the T-periodic solutions are expected to be close to the stable feasible fixed point, but 
the upper and lower bounds deviate from the stable fixed point by a factor of cT about the fixed point. As 
a result, we derive tighter bounds for y (t) in the following theorem. 

Theorem 5 Alternatively, y satisfies the bounds: 



ab-^cT aexp 



y (t) > max I ^-,5^, J - mm ^ 2^(,eT_i) ^ (T^^ W 

y it) < min I j^^^y !!!Kziz) j _ n,ax ^L^, (tS^) ) ' (4-9) 

where v^ax = max ^^'^^ G (t, s) f (s) ds, u„„ = mm G (i, s) / (s) ds, and G{t, s) = is the 

Green's function. 

Proof. The alternative bounds given by (|4.8p and (|4.9p are found working directly with the second equation 
of system (|2.4p . Given a T-periodic solution hit), it is known that the unique T-periodic solution of the 
first-order linear equation y' + cy = h(t) is 

rt+T 

y = J G{t,s)h{s)ds. 
Then, from the second equation of system (|2.4p , we have 

r-t+T 



y 



G{t,s) [dexp(2z) - f{s)]ds. 



Tight upper and lower bounds are derived for each of the two terms in the above integral considered sep- 
arately. By appropriately combining these bounds, we obtain (|4.8p and (14. 9p . See Appendix [El for the 
complete proof. □ 



The bounds in (14. 8|) and ()4.9p are tight for small forcing amplitudes. Combining these with the second 
inequality in ()4.7p . the bounds are tight for all forcing amplitudes. 
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Remark 2 Since y represents a squared relative velocity, a useful lower hound for y should he positive. In this 
sense, the second hound given hy (|4.8|) has the advantage that the lower hound is always positive if the external 
force f is not too large. In the model case f{t) = Aq cosiiot), it is easy to compute Wf'^Wi = ||/~||i = '^Aq/uj, 
II /II 2 = ^o-\/^, o,nd the hounds are direct. 

Having completed the mathematical analysis of the system (jl.2p . in the following section, we consider 
examples with realistic values of the parameters of the system representing different regions of the Saturn's 
ring and derive useful insights pertinent to the ring dynamics. 



5 Insights on Saturn's ring dynamics 

We begin with some corollaries of the theorems presented in the previous sections for the Saturn's ring 
system. 



Remark 3 The non-trivial, feasihle stahle fixed point for the unforced dynamic system is (Af, v'^^ 



' t)"^°"('M )t"' ' T^J • When externally forced hy periodic perturhations from the Saturn's moons, the mean 
aggregate mass and the velocity dispersion demonstrates periodic aggregation- disaggregation cycle with the 
same period as that of the moon's forcing. 

Corollary 1 (of Theorem\^ There exists a unique, asymptotically stahle limit-cycle for the (M ,v'^^^) system 
if the forcing period satisfies the condition 

1+4t2(1-£2)2 

where the right-hand- side depends on the mean optical depth and the normal coefficient of restitution corre- 
sponding to the specific location in the Saturn's ring system. 

The realistic parameters corresponding to the Saturn's ring system are e = 0.6, Tgyn = TLrb = T = 1. 
The optical depth is r = 0.1 for the F-ring and r = 1.5, 1, 0.5 for the B-ring outer edge. Hence, our result 
is applicable to the F-ring and also to the B-ring outer edge when r = 0.5. 

Now, we consider a systematic study of each region in the Saturn's rings, and start with the dynamics 
of the B-ring outer edge. 

Corollary 2 For a sinusoidal forcing function, f{t) — 2^cos(^^), the upper and lower bounds of the 
T -periodic solution for the mass of the aggregate is 

-Tg JUf (f\ 7,fS 

e7f(T^ < ''"'y-) _ <pv^(i-.^) f5Jl^ 

- 2v/4TT,to (1 - £2) - 

For v'^^-y (t), the hounds are obtained from Theorems^ and\^ with a = b = ^^^2 — 4t, Vcsc (Mq) — 0.5m/sec, 

Vth = Im/sec, Vmax = -Wrmn = ^^^T^ , II/II2 = vTLJ/WI^ = P , II/+II1 = /Lo 1/^ W I '^^ = 

2/3, , II/+II, - /f^o \f- {t)\ dt = 2/3, and \\f\\^ = 

We begin with studying the impact of varying the stirring rate /stir = on the amplitudes of the 

stir 

periodic solutions of (M (t) , ti^^,) (<)) , and Figures [Tal - [Tel plot the maxima and minima of the T-periodic 
solution of the (Af, u^^j) system obtained using numerical simulations against the derived upper and lower 
bounds for the T-periodic solution of the (Af , wj^^j) system from Theorem |4] and Theorem [5l The system 
parameters corresponding to the Saturn's ring system are e — 0.6, Tgyn = T^rh [T = 1), /3 — 0.5 and 
r = 1.5, 1, and 0.5 for Figure [Tal Figure ITbl and Figure [Tel respectively. The Saturn's B-ring including the 
Janus 2:1 density wave, and the B-ring outer edge have 1.5 and 1 as the typical values for the optical depth. 
Further, r = 0.5 is the typical optical depth value for the Saturn's A ring, Mimas 5:3 density wave and the 
Janus 6:5 density wave. 
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From (|5.1ip , it is clear that the amphtude upper and lower bounds for the periodic solutions of M [t) 
vary log- linearly with JsUr, and hence the corresponding curves in the plot are straight lines. It is interesting 
to note that the actual values for the maximum and minimum amplitudes (obtained through simulations) 
also appear to have a log-linear relationship with fsur with the same slope as that of the upper and lower 
bounds. As a result, an exponential increase in fsUr causes an exponential decay in the mean aggregate 
mass of the ring particles. This effect can be understood as follows. Large fsUr causes the rate of change of 
y (t) to be positive, and once v^^^ (t) > (same as y > ^ in (|2.4I) '). this causes a negative slope for M (t) , 
and essentially causing the exponential decay of M (t) . Further, the parameter fsUr has a negligible effect 
on the amplitudes of vf^^ (t) , as is evident from the negligible slope in the figure above. Notice that both 
the axes are in log-scale and the plots for the lower-bound of v'^^^ (t) does not appear in the plot because it 
takes negative values for all values of fsUr ■ 

We see that the effect of changing r only changes the magnitudes of maxima, minima, and the bounds 
on M [t) and v'^^-^{t). Otherwise, M (t) has a log-linear relationship with fsUr, with the same slope for 
all values of r. Comparing Figures [la] - [Tel the magnitudes of maxima M (t) increases with increase in 
r, and the minima decreases with increase in r. Further, the magnitude of the maxima of w^^j (t) remains 
unchanged with increase in r and the magnitude of the minima decreases with increase in r. The physical 
interpretation is that larger optical depth leads to a larger mean aggregate mass, but not large velocity 
dispersion. This is because the velocity dispersion is pegged near the fixed point value given by the strong 
erosion for v'^^y (t) > v^^^, which reverses the growth of M (t). For the same reason, increasing /stir reduces 
the upper value for M {t). 

Next, we study the impact of varying the forcing amplitude /3 on the periodic solutions of (M {t) , v^^^ (t)) , 
and Figures [2a| - [5c| plot the maxima and minima of the T-periodic solution of the (M, system obtained 
using numerical simulations against the derived upper and lower bounds for the T-periodic solution of the 
(M, w^pj) system from Theorem [4] and Theorem [5] The system parameters are the same as before and the 
stirring period is assumed to be equal to the orbital period of the ring particles around Saturn, in other 

words, /stir = 1- 

Notice that for small forcing amplitudes, the periodic solutions for M (t) and v'^^^ (t) has a small amplitude 
about their average values, which in this case, are the asymptotically stable fixed points (see Remark [3]) of 
the corresponding autonomous system of differential equations. With increase in /3, we see an increasing 
deviations in the amplitude of M (t) about the average values, while still non-negative, since the lower-bound 
from (|5.11|) is non-negative. Translating to the physical phenomena, the ring particle aggregates tend to 
achieve large masses as well as negligibly small masses, with increase in the forcing amplitude, in each period. 

The increase in forcing amplitudes causes a similar effect on "^^61 (^) M{t). But increasing the 

forcing amplitude beyond a certain threshold causes u^^j (t) to take negative values, which is non-physical. 
This sets a strict upper limit on the forcing amplitude, which in the case of the above figure, happens around 
P = 0.6. Comparing Figure [^aj and [5cl as r increases from 0.5 to 1.5, M (t) maxima increases, v'^^^ (t) maxima 
changes only slightly. 

Finally, notice that the upper and lower bounds for M (t) and v'^^^ (t) are tight for small forcing amplitudes. 
This means that the behavior of the system is well-characterized by the bound derived in Theorem [Jj and 
Theorem [5] for the regions in the Saturn's ring system where the infiuence of the moons are negligible. The 
analytical results provide a reasonable idea about the nature of the system and parallel the explicit numerical 
simulations. 

6 Conclusions 

We consider the mathematical model developed in [3] for understanding the dynamics of certain specific 
locations in the Saturn's ring system due to the periodic perturbations caused by the nearby moon. This 
model relates the mean aggregate mass and the velocity dispersion of the ring particles at a certain location 
in the Saturn's ring system in a manner similar to the popular predator-prey model in ecology, although it 
represents a totally different physical phenomenon. Further, this has been studied purely using numerical 
simulations so far. 

In this paper, we conduct a rigorous mathematical analysis of the dynamical system and derive results 
pertaining to the existence, uniqueness and stability of its solution. We analytically prove that the dynamical 
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stirring frequency relative to orbital frequency (f ^ in log-scale stirring frequency relative to orbital frequency (f ^ in log-scale 

(a) (b) 




Figure 1: Derived bounds compared to numerical simulations for varying /stir^ The system param- 
eters corresponds to the B-ring outer edge where the normal coefRcient of restitution e = 0.6, the forcing 
period of the moon is equal to the orbital period of the ring particles about Saturn (T = 1) and the forcing 
amplitude (3 = 0.5. (a) The mean optical depth is r = 1.5, (b) the mean optical depth is r = 1 and (c) the 
mean optical depth is r = 0.5. 
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Forcing ampUtude (P) in log-scale Poncing amplitude (p) in log-scale 

(a) (b) 




iV 10"' 10° 

Forcing amplitude (P) in log-scale 

(c) 

Figure 2: Derived bounds compared to numerical simulations for varying forcing amplitude 

/3: The system parameters corresponds to the B-ring outer edge where the normal coefficient of restitution 
£ — 0.6, the forcing period of the moon and the stirring period are both equal to the orbital period of the 
ring particles about Saturn (T = /stir = !)• (a) The mean optical depth is r = 1.5, (b) the mean optical 
depth is T = 1 and (c) the mean optical depth is r = 0.5. 
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system causes a periodic solution for the mean aggregate mass and the velocity dispersion of the ring particles 
at a period equal to the forcing period of the moon. Further, upper and lower bounds for the solution to 
the system are derived, that are especially tight for small perturbations caused by the moons. These bounds 
well-characterize the dependence of the system solution on various system parameters. Several useful insights 
are drawn about the Saturn's ring system based on the theoretical results presented here. 
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A Proof for Theorem [1] 

Deriving the first equation and inserting the second one, we obtain the following second order differential 
equation of Duffing type 

z" -\- cz' + 6dcxp(2z) ^ac + bf{t). (1.12) 

Now the periodic problems for (|1.2p and (|1.12p are completely equivalent. If (x, y) is a T-periodic solution 
of system (|1.2p with x{t) > for all t, then 2: = In a; is a T-periodic solution of (|1.12p . Conversely, if z is 
a T-periodic solution of eq. (|1.12p . then x ~ exp{z),y = (a — z')/b is a T-periodic solution of system (|1.2p 
with x{t) > for all t. 

Since ac > 0, eq. (|1.12p is of Landesman-Lazer type. Landcsman-Lazer conditions are classical and 
implies the existence of a periodic solution, see for instance |7i Theorem 1]. 

B Proof for Theorem [2] 

We know that system (|1.2p is equivalent to eq. ()1.12p , hence along this proof we will work directly with this 
last equation. 

Uniqueness. Assume that 21,22 are two T-periodic solutions of eq. (|1.12p . The difference S{t) — 
z\ (i) — Z2 (t) is a solution of a second order linear equation 

(5" c(5' ;i(i)<5 = 0, (2.13) 

where hit') = ^g; g'^p(2zi)-cxp(2z2) ^ consider two cases: 

• Case 1: 5{t) 7^ for every t: without loss of generality, we may take 5{t) > 0. By the Mean Value 
Theorem, 

h{t) = 26dexp(2^(i)), 
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where Z2(t) < £,{t) < zi{t) for every t. Hence, by using (|4.15l) . one obtains 



\\h\\^ = 2bd / cxp{2£_{t))dt < 2bd / exp{2zi{t))dt = 2acT. 
Jo Jo 

Using p.Sp . one may check easily that h e ^i.c, as defined in [6j Section 3]. Then, by [6j Corolllary 
2.5], the operator defined by the left-hand side of (|2.13|) is inversely positive, in particular (|2.13p is 
non-degenerate and its only periodic solution is the trivial one, which is a contradiction with the 
assumption that S(t) ^ for every t. 

• Case 2: 6{to) — for some to e [0,r]: by periodicity, there exists ti > to such that d{ti) — and 
5{t) 7^ for all t £]to,ti[. Again, we may assume without losing generality that S{t) > for all 
t G]to,ti[. Define the truncated function h{t) as h{t) if t G [to,^i] and if t G [0, io[U]ii, T]. By 
reasoning as before, h G ili^c- Then, as a consequence of the results in [B] (see Remark 2.2 therein), 
the distance between two consecutive zeroes of a solution of S" + c5' + h{t)5 = is greater than T, but 
(5 is a solution vanishing at toi^ii which is a contradiction. 

After this analysis, the only remaining possibility is d{t) — for every t, so the T-periodic solution is unique. 

Stability. The proof for asymptotic stability is also based on the results from [6j. Let z{t) be the unique 
T-periodic solution of eq. (jl.l2p . A standard computation gives that the topological index of z is 7(2;) = 1. 
The variational equation is u" -\- cu! -\- h(t)u = with h(t) = 2bdexp{2z(t)) and again h € ili.c as a result of 
p.5|) . Then the asymptotic stability follows directly from [B] Theorem 3.1]. 



C Proof for Theorem [3] 

We use the terminology and notation from [3]. A lower solution is found by fixing a constant a such that 

Mexp(2a) =ac + & 11/11^. (3.14) 

An upper solution is easily obtained as follows. Take F the unique T-periodic solution of z" + cz' = f{t) 
with mean value zero. Then, /3(t) = —M + F{t) with M big enough is an upper solution with /3 < a . 
To apply 131 Theorem 1.2], we have to check two hypotheses. The fist one is the finiteness of the number 
of T-periodic solutions, which is direct from the analiticity of the field on the state variables. The second 
hypothesis is condition (1.6) in 4 . Fix g{t,z) = bdexp{2z) — ac~ bf{t). For s e [/3(t),a], 

— g(i,s) = 26dexp(2s) < 26dexp(2a) = 2ac+2fe||/||^ < T + ^ 
by (|3.6p . which is just (1.6) in [4l. The result is done. 



D Proof for Theorem [4] 

We work again over the equivalent eq. (|1.12p . Integrating ()1.12p over a whole period 

bd [ exp{2z)dt = acT. (4.15) 

"'0 

By the integral mean value theorem, there exists to G [0, T] such that 

exp(2z(io)) = ^. (4.16) 

Observe that z (i ) — z (to ) G Hq (tp , T -|- tg ) , then by the Sobolev inequality with the optimal constant (see for 
instance [S]), we get 



\zit)-zito)\<^\\z'\\,. (4.17) 
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On the other hand, multiplying (jl.l2p by z' and integrating over [0,T] and applying Cauchy-Schwarz in- 
equality, we get 

c\W\\l = b f fit)z'dt<b\\.f\\Jz'\\^, 
Jo 

hence, 

Ik'll2<^ 11/112- (4.18) 

Now, combining (j4A6| . (|4Tf|) and (|4J8l) . 

and 

= zito) + z{t) zito) > 1 In (g) - . 

The first inequality of (|4.7p arises from here taking z = Ina; and finding the value of x. 

The bound for y is found as follows. Let such that = mintg[g j^j z(t). Integrating (jl.l2p on 



z'{t) + cz(t)-cz{t.,) + bd J exp{2z)dt^J ac + bf{s)ds < acT + b\\f+\\^ . 
Since cz{t) — cz{t^.) + bd exp{z)dt > 0, then 

z'{t)< ac + bf{s)ds < acT + b\\f+\\^. (4.19) 
Analogously, Let t* such that z{t*) — max^^^Q x] z{t). Integrating (|1.12p on [t* ,t], 

z'{t)+cz{t)~cz{t*) + bd f exp{2z)dt= f ac + bf{s)ds > -b\\f~\\^. 



t' 

Then, 

ft'+T 



pt pt +1 

z'{t)>~bdj exp{z)dt- b\\f^\\^> -bd J exp{z)dt - b\\f~\\^ . 

Since z{t) is T-periodic, /j*, ^'^ exp{2z)dt = exp(2z)dt, then by using (|4.15p . 

z'{t)>-acT-b\\f-\\^. (4.20) 
Finally, inserting (|4.19p — (|4.20|) into y — -^(0 — z'), one obtains the second inequality of (|4.7p . 



E Proof for Theorem [5] 

Note that 



<G(M)< '"^^'^^ 



exp(cT) - 1 - ^ ' ' - exp(cr) - 1 
for all t,se [o,r]. 

In the following steps, we derive upper and lower bounds for both the terms in the above expression, 
which will be used later to obtain bounds for y (t) . 
Bounds for u (t) — J^^^ G (t, s) exp (2z (s)) ds: 
Using the bounds for G {t, s), 

acT , , acTe^^ 

<^^t)<———. (5.21) 



bd (e^T - 1) - ' - bd {e^T _ i) ■ 
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The next set of bounds are derived by using the bounds derived for z (t) in the proof of Theorem U] and 
noting that J^^^ G {t, s)ds = ^ 

-exp y j < K(t) < -exp ( — ) . (5.22) 

Combining (IOT|) and (j??^ , we get 



( acT a ( 5Vr||/L\\ . / acTe"'^ a (b^/T\\f\L\\ ,^ 

[mW^Y m [ 7^ j j ^ " W ^ [mW^Y m [—V^ ) ) ■ (^-2^) 

Bounds for v (t) = /j+f G {t, s) f (s) ds: 

Let pit,s) = jt+?'^^'^^\^^^ ■ We know, f^^^ G{t,s)ds = i. Hence, v (t) = f^^f p (t, s) f {s) ds. Since 
fmin < Js=J p{t, s) f (s) ds < fmax, wc havc the first set of bounds 

fn^ ^ frnax_ (5.24) 

c c 



Similarly, using the Cauchy-Schwartz inequality, we have p (t, s) / (s) ds < II/II2 y /^^^ s) ds and 

hence 



11/11 J^^^^. (5.25) 



Another set of bounds based on the following inequality 



are 



t+T rt+T 

G{t,s) f - {s)ds <v{t) < / G{t,s)f+{s)ds 



—LLJW < (t) < (5.26) 

c c 

"^""1'''^ <t;(i) < M^JJi^. (5.27) 



And finally, if the v {t) is computable for all t, then, 

Vniin <v{t)< Vjnax- (5.28) 

Note that (|5.24p and (|5.26|) are equivalent. Hence, combining (j5.25p - (|5.28p . we get 

-iiriL -iiriiie-^ 

' (e^T-l) 



V (t) > min v> 



V it) < mm 11/11, J ^^^^^r_^y ' ^i^r^iy ' (5-29) 



Using and dOH), (SH) and gH) follow directly. 
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